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Skyrmions are the magnetic defects in an ultrathin mag- 
netic film, similar to the bubble domains in the thicker films. 
Even weak uniaxial anisotropy determines its radius unam- 
biguously. We consider the dynamics of the skyrmion decay. 
We show that the discreteness of the lattice in an isotropic 
2D magnet leads to a slow rotation of the local magnetization 
in the skyrmion and, provided a small dissipation, to decay 
of the skyrmion. The radius of such a skyrmion as a func- 
tion of time is calculated. We prove that uniaxial anisotropy 
stabilizes the skyrmion and study the relaxation process. 

PACS numbers 75.50. Ss, 75.70.Ak, 75.70.Kw, 75.40.Gb 



Skyrmions are topological excitations of 2D magnet or 
ultrathin magnetic films similar to the well-known bubble 
domains in the industrial magnetic memory materials |Q . 
The former play an important role as letters in magnetic 
records. A natural question arises whether the skyrmions 
can be utilized in a similar way. To answer this question 
one must adjust the existing theory to a realistic condi- 
tions of real films with their anisotropy, defects, dissipa- 
tion, and discreteness. It is necessary also to consider 
how the skyrmion can be created and destroyed, i.e. the 
skyrmion dynamics, ignored in the previous studies. In 
this paper we address these problems. 

The Skyrmion was first discovered by Skyrm pi who 
considered it as a localized solution in a model of the nu- 
clear matter. Belavin and Polyakov (BP) || have shown 
that the skyrmion is a topologically nontrivial minimum 
of energy for the so-called n-field, the classical continuous 
limit of the Heisenberg model. It realizes the mapping 
of the plane in which the spins are placed onto sphere of 
the order parameter with the degree of mapping 1. 

Though the n-field model was inspired by the stud- 
ies of magnetic films, the first more or less direct obser- 
vation of the skyrmions was made on the 2D electron 
layers under the QHE conditions B following an ear- 
lier theoretical predictions ||. An indirect observation 
of skyrmion effects in quasi-2D magnets was reported by 
F.Waldner Q who found the skyrmion energy from the 
heat capacity measurements in a good agreement with 
the theoretical prediction by BP. This experimental ob- 
servation is highly nontrivial since the real magnets are at 
least weakly anisotropic. Due to the existence of Gold- 
stone modes in the Heisenberg magnet even very weak 
anisotropy can change the excitations crucially. 



Let us first approach the problem with simple dimen- 
sionality arguments. The skyrmion is a static excitation 
of the homogeneous ferromagnetic state localized in a 
circle of the radius R. From the dimensionality consider- 
ation and from the BP results the skyrmion energy does 
not depend on its size and is equal to 47rJ|m| where m 
is the degree of mapping. Hence, the exchange energy 
1/2/ J(VS) 2 d 2 a: of a skyr mion does not depend on its 
radius R. The anisotropy energy 1/2/ XS^d 2 x of the 
skyrmion is proportional to R 2 and decreases together 
with R. Therefore the exchange energy of the fourth or- 
der in the derivatives 1/2/ «;(AS) 2 d 2 x is crucial for the 
skyrmion. It stabilizes the skyrmion if k > 0. Now the 
total energy of the skyrmion depends on its radius and 
has the minimum at R ~ (k/A) 1 / 4 . One can estimate n 
by the order of magnitude as Ja 2 where a is the lattice 
constant. Then R ~ (al) 1 / 2 <C l\ where l\ = \/J/X is 
the domain wall width. The energy of such a skyrmion 
differs by about XR 2 ~ \/kX from the classical skyrmion 
energy 4nJ 3> VkX. 

Let us consider the situation more closely. The classi- 
cal two-dimensional Heisenberg exchange ferromagnet in 
continuous approximation is described by the Hamilto- 



Ho = \ J J(VS) 2 d 2 x 



(1) 



with the constraint on the vector field S(r): S 2 (r) = 1. 
An obvious minimum of such a Hamiltonian is the ho- 
mogeneous ferromagnetic configuration in which all the 
spins are parallel S(r) = const. The simplest topologi- 
cally non-trivial minimum of the Hamiltonian (Q) is given 
by the Skyrm solution: 
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It describes a skyrmion of radius R with the center placed 
in the origin. The observation point is indicated by 
the polar coordinates r and <fi; tj) is an arbitrary an- 
gle. We have already mentioned that the energy E of 
the skyrmion (||) does not depend on its radius R It also 
does not depend on the angle ip. 
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Any static distribution of magnetization S(r) satisfies 
the equilibrium equation: 
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(3) 



The r.-h.-s. of eqn. (|3|) is added as a Lagrangian factor 
which ensures that S 2 (r) = 1 at any point r. For H = 
H() + Hi where Hi is a perturbation we look for a solution 
in a form S = So + Si, where So(r) is determined by eqn 
(Q) and Si(r) is perpendicular to So and satisfies the 
linearized inhomogeneous equation: 
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s=s 



where the tensor kernel Kij(r,r') of the linear operator 
K is given by 

K ij (r,r') = -JS{r-r') x 

(5ijA r i - Sij (S • AS ) - S iS 0j A r i). 

The derivatives dSo/dR and dSo/dip are the zero 
modes and, hence, satisfy the homogeneous equation 
KdSo/dR = KdS /dip = 0. Therefore the right hand 
side of eqn. (0) must be orthogonal to the vectors 
dSo/dR and dSo/dtp. The two orthogonality condi- 
tions allow us to determine both R and ip fixed by the 
small perturbation Hi. Considering a special perturba- 
tion Hamiltonian 



Hi = 



\ J ( K (AS) 2 + A(1-S 2 )) d 2 x, 



we find from the orthogonality condition: 
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(5) 
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where L = log[(3A/8K) 1 / 4 ^]. Logarithm in eqn. @ 
comes from the divergent integral / (1 — Sq z )8Sq z / dRd 2, x. 
It was cut off at the radius r — l\ at which the perturba- 
tion theory fails. Due to the axial symmetry ip remains 
zero mode. Note that the skyrmion does not exist for 
k < 0. There is no general argument in favor of positive 
k, but for standard discrete Heisenberg model n is posi- 
tive. Now we are in position to consider the dynamics of 
the skyrmion. 

The dynamics of the unit vector field S(r, t) is given 
by the Landau-Lifshitz equation R]: 



S(r,t) = -gS(r,t) x 



vS(r,t) x 
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It can be checked that at v = the equation of motion 
(0) conserves the magnetization and energy of the field 



S(r) as well as the local constraint S 2 (r) = 1. A small 
dissipation term at v ^ allows for the relaxation pro- 
cesses. 

We consider only the slow skyrmion dynamics. It 
means that we present again the Hamiltonian H as a 
sum Ho + Hi, where Hi is a small perturbation to the 
exchange Hamiltonian Ho . We are looking for a solution 
in the form S = S + Si, S • S x = 0, where |Si| « |S | 
and So(r; R(t),ip(t)) is the standard skyrmion solution 
(^) with the parameters R and ip slowly varying in time. 
We also consider the dissipation and the term S as a 
perturbation. Substituting only So(r; R(t), ip(t)) in the 
perturbation terms in eqn. (R) and requiring the orthog- 
onality of the perturbation terms to both zero modes 
8S Q /dR and dS /dip, we obtain equations of motion for 
R and ip: 
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(8) 
(9) 



where uj = ip and R is a scale at which the perturbation 
expansion breaks down, i.e. |So(r = R)\ ~ |Si(r = R)\. 
Eqns (^) and (^|) have a fixed point with ui — and 
R = Ro, where Ro is determined by eqn. (|]). 

If R deviates slightly from Rq so that AR = R — Rq 
is still small enough, one can put R = l\ into (^|) and 
obtain w = AXgAR/R. Using (|), we find: 



AR = -AXvAR. 



(10) 



"Small enough" AR means that we still can use R~l\. 
However, if AR is large, the cut-off scale is determined 
by the finite frequency R~ 1^ — y/gj/\uj\. Eqn. ( |To| ) is 
valid if h < L or AR/R < 1/4. 

In the opposite case l u -C l\ and R ^> Ro one can 
neglect the second term in (^) . Then equations of motion 
(||, 0) read as follows 



= Xg: R = —vXR. 



(11) 



For l u <c l\ and R <C -Ro one can drop the last term in 
(0). Then: 
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(12) 



Thus, the easy-axis anisotropy together with the 
fourth-order exchange term fix the radius (||) of the 
skyrmion if k > 0. In the opposite case K < there is 
no stable configuration unless the higher order exchange 
interaction is taken into account. We also have shown 
that the skyrmion reaches its equilibrium radius with the 
characteristic time t r = (vX)~ 1 . 

Another property of a real film which should be taken 
into account is the lattice discreteness. We have men- 
tioned earlier that in the continuous model the skyrmion 



2 



is a topological excitation and as such cannot dissipate. 
However, in the discrete lattice the continuity of the field 
S(r) is lost and the very notion of the topological ex- 
citation becomes inconsistent. Therefore, the skyrmion 
configuration in the discrete lattice is unstable. Moreover 
it is sufficient to remove one plaquette in the center of 
the skyrmion to make it unstable, (see for example 0). 
We will imitate the discreteness effect by considering a 
hole in the center of the skyrmion. In the picture due to 
BP the skyrmion is described by a meromorphic function 
which has a pole in the center of the skyrmion. This pole 
cannot be rid off by any continuous change of the field 
S. To allow a skyrmion to dissipate we should punch a 
small hole in its center. 

Let us first consider the skyrmion without anisotropy. 
"Punching a small hole in the center" means the substi- 
tution J — * J8(r — r ) in cqn. (|l|) where r <C R is the 
radius of the hole and 6(x) is the step function: 9(x) = 
for x < and 0{x) = 1 for x > 0. Hence, the perturba- 
tion Hi is given by: 



H, = -J 



(r -r)(VS) 2 



(13) 



Employing the orthogonality condition of the right hand 
side of the linearized eqn. (fj]) to the zero modes and 
neglecting all the terms of the higher order in ro, one 
finds that eqn. (||) still holds, but eqn. @ must be 
replaced by: 
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where R « l u = ^JgJ/\w\ is the scale where the per- 
turbation scheme breaks down. With the logarithmic 
accuracy one can write ln(i? 2 /R 2 ) — ln(i? 2 /r 2 ,). Finally, 
substituting to from (O) into (pi), we find: 



R 3 R = 



In 



(15) 



From this equation one can conclude that the skyrmion's 
life-time is roughly proportional to its radius in the fourth 
power. 

Returning to the field with the easy-axis anisotropy 
and the fourth-order term, let us introduce again a hole 
in the center. In this situation the perturbation is given 
by the sum of the two terms (||) and (|IJ). Using the 
same orthogonality trick one gets: 



(0) and ([TJ) are valid as well after the same substitu- 
tion. In the case k < (3/8)Jr-Q, however, the effective k 
is negative and the stable skyrmion exists no longer. 

In conclusion, we have shown that in a real ultra- 
thin ferromagnetic film with easy-axis anisotropy the 
skyrmion known for the isotropic model still exists, but it 
acquires a definite radius Ro given by eqn. (||). By the or- 
der of magnitude Rq ~ \Jal\ <C I \ where I \ is the domain 
wall width. Once we made a domain with the reversed 
magnetization in a ferromagnetic film it shrinks down to 
the size Rq ~ Inm. The discreteness of the lattice in the 
isotropic model leads to a finite skyrmion life-time which 
is roughly proportional to the fourth power of its radius. 
However, anisotropy together with the higher order ex- 
change interaction stabilizes the skyrmion. At finite tem- 
perature it can decay through an instanton configuration. 
Our results allow to understand why the activation en- 
ergy found by Waldner is so close to AttJ: the difference 
is expected to be of the relative order a/l\ ~ 10~ 2 . The 
detailed thermodynamics of skyrmions will be published 
elsewhere. 
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After substitution k — > k = k — (3/8) Jr^ eqn. (|FoJ) ac- 
quires the same form as eqn. ([}]). It means that as long 
as k > (3/8)J7q the skyrmion is stable and its radius is 
defined by eqn. (ph with k instead of k. Equations (10), 
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